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Abstract 

In M. d , for any dimension d > 1, expansions of self-intersection local 
times of fractional Brownian motions with arbitrary Hurst coefficients 
in (0, 1) are presented. The expansions are in terms of Wick powers 
of white noises (corresponding to multiple Wiener integrals), being 
well-defined in the sense of generalized white noise functionals. 
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1 Introduction 



Intersection local times of fractional Brownian motion (fBm) has been studied 
by many authors, see e.g. the works done by Gradinaru et al. [GRV03] . 
Nualart et al. |HN07j , [HN05] , Rosen [Ros87] . and the references therein. 

We may consider intersections of simple paths of a fBm with themselves 
or with other independent fractional Brownian motion as in [NOL07j . Within 
the white noise analysis framework, we concentrate here in the simple self- 
intersection problem using a similar approach to Faria et al. [dFHS W97j . 
[dFDSOO] to study the case if = 1/2 (the classical Brownian motion). This 
approach has the advantage that the underlying probability space does not 
depend on the Hurst coefficient under consideration. As a consequence, we 
analyze the self-intersections of a fBm, with Hurst coefficient being any possi- 
ble value in (0, 1). As one may expect, independently of the Hurst coefficient 
under consideration, self-intersections become scarce as the dimension d in- 
creases. 

An informal but suggestive definition of self-intersection local time of a 
generic d-dimensional fBm B is in terms of an integral over Donsker's 5- 
function 



intended to measure the amount of "time" the process spends intersecting 



A rigorous definition, such as, e.g., through a sequence of Gaussians ap- 
proximating the 5-function, 



will make L increasingly singular, and various "renormalizations" have to be 
done as the dimension d increases. For d > 1 and l/d<H<3/(2 + d), the 
expectation diverges in the limit and must be subtracted. The L 2 -properties 
related to this situation have been analyzed in |HN05j . For other H values 
further kernel terms must be subtracted (Theorem [7| below). 

In this work we are particularly interested in the chaos decomposition 
of L. We expand L in terms of Wick powers [HKPS93J of white noise, an 
expansion which corresponds to that in terms of multiple Wiener integrals 
when one considers the Wiener process as the fundamental random variable. 
This allows us to derive the kernels for L. Due to the local structure of the 



L 




itself. 




e > 
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Wick powers, the kernel functions are relatively simple and exhibit clearly the 
dimension dependence singularities of L (Proposition [S]). For comparison, we 
also calculate the regularized kernel functions corresponding to the Gaussian 
5-sequence mentioned above. 

The paper is organized as follows. In Section 2 we review the necessary 
background of white noise analysis and in Section 3 we present the main re- 
sults of this work and their proofs. We shall mention that after the conclusion 
of our work we found the recent paper [Rez07] with related results. 

2 Gaussian white noise calculus 

In this section we briefly recall the concepts and results of white noise analysis 
used throughout this work (for a detailed explanation see e.g. |BK88j . [Hid75] . 
[HKPS93j . [H0UZ96] , [Kon75j . [KT80aj . [KT80bj . [Kuo96j . [Oba94j ). 

2.1 Fractional Brownian motion 

The starting point of white noise analysis is the real Gelfand triple 



where L 2 d {R) := L 2 (R, R d ), d > 1, is the real Hilbert space of all vector valued 
square integrable functions with respect to the Lebesgue measure on R, S d (R) 
and S' d (R) are the Schwartz spaces of the vector valued test functions and 
tempered distributions, respectively. We shall denote the L^(M)-norm by | • | 
and the dual pairing between S' d (R) and S d (R) by (-,-), which is defined as 
the bilinear extension of the inner product on L d (R), i.e., 



for all g = fox, ...,&,) G L 2 (R) and all f = (A,. ..,/,) G S d (R). By the 



Minlos theorem, there is a unique probability measure /x on the a-algebra B 
generated by the cylinder sets on S^(R) with characteristic function given by 



S d (R) C L 2 d (R) c S' d (R), 





feS d (R). 
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In this way we have defined the white noise measure space (S' d (JNL),B, fi). 
Within this formalism, a version of the <i-dimensional Wiener Brownian mo- 
tion is given by 

B(t) := (fa, l m ), (co d , l m )) , (u u uJ d ) G S' d (R) 

where 1^ denotes the indicator function of a set A. For an arbitrary Hurst 
parameter 0<H<1, H^~, a, version of a ^-dimensional fractional 
Brownian motion is given by 

B H (t) := ((u 1 ,M H l [0!t] ),...,(u; d ,M H l m )) , K,...,^) G S' d (R), 

where, for a generic real valued function /, and for ~ < H < 1, 

(M H f)(x) := n / dyf(y)(y-x) H -l, (1) 



r(H-i) 

provided the integral exists for all iGl, while for < H < ^, 

| - H)K„ /•» /|» - /(! + y) 



provided the limit exists for almost all x G R (for more details see e.g. |Ben03] 
and [PTOO] and the references therein). Independently of the case under 
consideration, the normalizing constant Kh is given by 

k »^(h + \) (4 + p s ( (1 + s) ^- s ^))-\ 

Apart from the operators ([I]), ([2]), we shall also consider the operator 
defined for | < < 1 by 

(M+f)(x) := * H [ X dyf(y)(x - y)*~\, (3) 

provided the integral exists for all x G R, and the operator defined for < 
H < \ by 
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provided the limit exists for almost all 

There are several examples of functions / for which M^f and M^f exist 
for any H G (0, 1). For instance, / = l[o,t], t > 0, or / G Si(R). For functions 
fi, /2 being either one of these two types it is easy to prove the following 
equality 

[ dxh{x){M H f 2 ){x)= [ dx{M+h){x)h{x), 
Jr Jr 

showing that Mh and are dual operators. For more details and proofs 
see e.g. [Ben03j and the references therein. 

2.2 Hida distributions and characterization results 

Let us now consider the complex Hilbert space L 2 (S' d (M), B, fi). This space 
is canonically isomorphic to the symmetric Fock space of symmetric square 
integrable functions, 

L 2 (S' d (R),B,fi) ~ (0SymL 2 (R fe ,A;W fc x)j , 
leading to the chaos expansion of the elements in L 2 (S' d (R), B, fi), 

F(U U W„) = ^ U ? H1 : ® • • • ® : "T* f{nu..,n d) ), (5) 

(ni,...,n d )eN d 

with kernel functions f( ni ,...,n d ) i n the Fock space. For simplicity, in the sequel 
we shall use the notation 

d d 
n = (m, • ■ ■ ,n d ) G N d , n = ^2rii, n! = JJni!, 

8=1 8=1 

which reduces expansion (jSJ) to 

F( W )=^(:^:,f n }, uES' d (R). 

n€N d 

To proceed further we have to consider a Gelfand triple around the 
space L 2 (S' d {R), B, //). We will use the space (S)* of Hida distributions 
(or generalized Brownian functionals) and the corresponding Gelfand triple 
(S) C L 2 (S' d (R),B,fi) C (S)*. Here (S) is the space of white noise test 
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functions such that its dual space (with respect to L 2 (S' d {R), B, //)) is the 
space (S)*. Instead of reproducing the explicit construction of (S)* (see 
e.g. [HKPS93J), in Theorem [2] below we define this space by its S'-transform. 
We recall that given a f G S'd(R), and the Wick exponential 

: exp«u;,f» : := £ I<: cu® n :,f°») = C{f)e™ , 

n€N d 

we define the S'-transform of a $ G (S)* by 

5$(f) :=«$, :exp((-,f)):)), Vf G ^(R). (6) 

Here ((•,•)) denotes the dual pairing between (S)* and (S) which is defined as 
the bilinear extension of the sesquilinear inner product on L 2 (S' d (M), B, [/,). 
We observe that the multilinear expansion of (jSJ), 

5$(f) :=£(F n ,f® n >, 

n 

extends the chaos expansion to $ G (S)* with distribution valued kernels F n 
such that 

(($,^>> = $>!(F n ,y, n >, 
n 

for every generalized test function (p G (S) with kernel functions ip n . 

In order to characterize the space (S)* through its ^-transform we need 
the following definition. 

Definition 1 A function F : S^R) — > C is called a U -functional whenever 
1. for every f 1; f 2 G S'd(R) the mapping R 9 A i — > F(Afi + f^) /ias an entire 
extension to A G C, 

U. iaere are constants K\,Ki > swca iaai 

\F(zf)\ < Kiexp (K 2 \z\ 2 1| ff) , G C,f G S^R) 

/or some continuous norm ||-|| on Sd(R). 

We are now ready to state the aforementioned characterization result. 

Theorem 2 ( [KLP + 96j . |PS91j ) The S -transform defines a bijection between 
the space (S)* and the space of U -functionals. 
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As a consequence of Theorem [2] one may derive the next two statements. 
The first one concerns the convergence of sequences of Hida distributions 
and the second one the Bochner integration of families of the same type 
of distributions (for more details and proofs see e.g. |HKPS93j . |KLP+96j . 



Corollary 3 Let ($ n ) ne N be a sequence in (S)* such that 

(i) For allieS d (R), ((5$ n )(f)) neN is a Cauchy sequence in C, 

(ii) There are K\,Ki > such that for some continuous norm ||-|| on Sa(R) 



Then ($ n )„ g N converges strongly in (S)* to an unique Hida distribution. 

Corollary 4 Let (Q, B, m) be a measure space and A t— > $a be a mapping 
from Q to (S)*. We assume that the S -transform of $a fulfills the following 
two properties: 

(i) The mapping X \— > (S$\)(f) is measurable for every f G Sd(K), 

(ii) The S&\ obeys a U -estimate 



for some continuous norm \\ ■ \\ on Sd(R) and for C\ G L 1 (fi,m) ; C<i G 



[PS93]). 



one has 



\(S$ n )(zf)\ < K ie K2 ^ m \ Vz G C,f G S d (R),n G N. 



\(S$ x )(zf)\ <C x {\)e 



C 2 (A)| 2 | 2 ||f| 



zeC,i e S d (R) 



L°°(fi,m). 



Then 




and 



S 



( [ dm{\) $ A ) (f) = / dm{\) 5$ A (f) 



\Jn J Jn 
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3 Intersection local times 

Proposition 5 For t ^ s the Bochner integral 

5(B H (t) - Bff(a)) := ( i-V / dA e «CBH(*)-BH(.)) (7) 

\2ttJ J Rd 

is a Hida distribution with S-transform given by 



S5(Bff(t)-Btf(s))(f)=(-^ -) e ai* — i 



'27T t - S 



.u_,|2H 



2 



(8) 

/or all f G Sd(R) . 

Proof. The proof of this result follows from the application of Corollary H] to 
the S'-transform of the integrand function in ([7|) with respect to the Lebesgue 
measure on M d . 
Since 

Se a(B H (t)~V H (s))tf\ = e -2f(t-s) 2H e i\f R dxf(x)(M H l [Sit] )(x)^ ,g\ 

for t > s, and 

Se i\(B H (t)-B H (s))^ _ e -l>f( s - t ) 2H e -iXf M dxt(x){M H l [tA )(x)^ ^ 

for t < s cf. e.g. [H KPS93] . clearly in both situations the measurability 
condition is fulfilled. 

Independently of the particular case under consideration, we observe that 
for all z G C and all f G S^fKO we find 



|£ e iA(B H (t)-B ff ( S ))^ f) | = 



Al 2 

- J-J- \t-s \ 2H ±iz\ J m dx f (x) (Mb 1 [ S At ,s vt] ) O) 



^|t-s| 2H p -^|t-^| 2H + l^l|A||/ M ^fW(A/ ff l [sAt , flVt] )(x)| 



< e" T - |C-S| e 

where the latter exponential is bounded by 

i i 2 i 

g |f_J|2H | /k da: f ( x ) ( M » 1 [« At ,s vt] ) (a>) 
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because 



IAI 



t - s\ 2H + \z\\\ 





z 




\t - 


- s 


H 



— It - s 
2 1 



t-s 



2H 



As a result, 



dxf(x)(M H l [sAtiSVt] )(x 
dxf(x)(M H l[a/\t,svt])(x 
dxf(x)(M H Ti-[sAt,svt])(x) 



H 



|^ e iA(B H (t)-B H ( S ))^ f ^| < e -^| 4 _ s |2H e ^7i27/|/H da;f ( a; )( M ff 1 bA t , s v t ])W| 

where, as a function of A, the first factor is integrable on R d and the second 
factor is constant. 

An application of the result mentioned above completes the proof. In 
particular, it yields (jHJ) by integrating (jHJ), (TTO]) over A. ■ 

In order to proceed further the next result shows to be very useful. It 
improves the estimate obtained in }Ben03t Theorem 2.3] towards the charac- 
terization results stated in Corollaries [3] and HI 

Lemma 6 Let H 6 (0, 1) and f G Si(K) be given. There is a non-negative 
constant Ch independent of f such that 



< C H \t - s\ fsup |/(x)| + sup + 

\a;GM xSK 



dxf{x){M H \ [sA ){x) 
for all s < t. 

Proof. Since Mh and are dual operators, 

/ dxf(x)(M H l M )(x)= f dx(M+f)(x), 
Jm. J s 

and thus 



dxf{x){M H %,t]){x) 



< |t- S |sup|(M+/)(x)|. 
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It has been shown in |Ben03t Proof of Theorem 2.3] that for < H < ~, 

sup \(M+f)(x)\ < C' H f^ 2 — sup |/(x)| + ~77~~T sup \f'(x) 

for some constant C' H > independent of /, while for ~ < H < 1, 
l(M+/)(x)| 
< C£ / rfy |/(y)||j, -afH 



= C£/ dy|/(y)||y-x| H -«+C^/ dy|/(y)||i/ 

^|j/-x|<l ^|j/-x|>1 

< SUP + C'k I dy\f{y)\\y-x\ H -l 

tl — 2 ./|j/--x|>l 

for some constant > also independent of /. Concerning the latter 
integral, observe that 

1/2 

dy\f(y)\\y-x\ H -3 < \f\ ( / dy\y x\ 21 



y— x\>l \J\y—x\>l 

' 1/1, 



leading to 



SUp | (M+f) (*) | < C'k f-^-y SUp + - ? =i=|/| > ) • 



We are now ready to state the main result on intersection local times Lh 
as well as on their subtracted counterparts For simplicity we shall use 

the notation 

A := {(t u t 2 ) G M 2 : < h < t 2 < !}■ 

Theorem 7 For any H e (0, 1) and for any pair of integer numbers d > 1 
and N > suc/i #icrf 2N(H — 1) + rfiJ < 1, the Bochner integral 



LyP := / d 2 t5^(B H (t 2 )~B H ( tl )) 
J A 



is a Hida distribution. 
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Proof. To prove this result we shall again use Corollary H] with respect to 
the Lebesgue measure on A. For this purpose let us denote the truncated 
exponential series by 



exp N {x) 



EX 



n=N 



It follows from (jBJ) that the ^-transform of 5^ is given by 
5(«jW(BH(t 2 )-B H (*0))(f) 



(11) 



1 



/ 27r|t 2 - h 



H 



exp N 



2\t 2 -tx 



2H 



dxi(x)(M H \ [tlM ){x) 



which is a measurable function. 

In order to check the boundedness condition, on Sd(R) let us consider the 
norm || • || defined by 



|f|| := [J2 (mp\f i (x)\+ sup 1/^)1 + 1/, 



2\ 2 



f=(h,...,f d )eS 6 



(12) 



We observe that for d — 1 this norm reduces to a continuous norm on Si ( 
= sup|/(x)|+sup|/'(x)| + |/|, feS 1 (R), 



which implies the continuity of the norm (fT2"|) for higher dimensions. By 
Lemma Owe obtain 



dx f(x)(M H l[ tut2 ])(x) 



211*112 



<Cji\t 2 -h\ z \\f 



and thus 



< 



1 V (C 2 



H I |2u i 1 2— 2H||*||2 

, , ... , ex P7V — \z\ Ih-hl ||f" 

27T t 2 — *1 1 I V 2 



for all z G C and all f e Sd(R). Estimating the function exp^ by 



expjv 



r<2 



_ , |2iV(l-//) p ^| Z p||f||2 
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we then obtain 

d 



\S(5^(B H (t 2 )-B H ( tl )))(zf)\ < (J=\ \t 2 - tl \ 2N ^- dH e^ 2 ^\ 

where \t 2 -t x \ 2N ^~ H) ~ m is integrable on A if, and only if, 2N(l-H)-dH > 
— 1 . The proof is completed by an application of Corollary HI ■ 

As a consequence, one may derive the chaos expansion for the (truncated) 
local times 

Proposition 8 Given a H G (0, 1) and a pair of integer numbers d > 1 and 
N > such that 2N(H — 1) + dH < 1, the kernel functions Fjj n are given 
by 

1 ( 1 \ d ( l\ n f 1 2n 

F H , 2n {u X) ... : u 2n ) = -{-j=^ (--J JjH lh _ tiimn+dH U( M ^M)( 

for each n G N d such that n > N. All other kernel functions -F# n are 
identically equal to zero. 

Proof. By Corollary HI the S'-transform of the (truncated) local time is 
obtained by integrating (jTTj) over A. Hence, given a f = (f\, fd) G Sd{ 



N 2 n \t 2 -t x \ 2Hn 

n=N ' 1 



X 



E wl .!. wi n( / dx f,(x)(M H l [tlM )(x)) 

a>£(-j) 



2/i 



1 

\2Hn+dH X 



n=JV 

d / „ \ 2n, 



I ,"■ * — 1 \v R 



n\,— ,n d j = i 
n±-\ \-7i d — n 



Comparing with the general form of the chaos expansion 

4f = £<:^:,i^ n >, 
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one concludes that 

1 / 1 \ d f l\ n f 1 2n 

F HM (u 1 ,...,u 2n ) = -\^=j JjH ^_ t ^ Hn+dH Y[(M H l ltlM )( Ui 

for each n G N d such that n > N, while all other kernels Fh,h vanish. ■ 

Theorem [7] shows that for d — 1 all intersection local times Lfj are well- 
defined for all possible Hurst parameters H in (0, 1). For d > 2, intersection 
local times are well-defined only for H < 1/d. Informally speaking, for H > 
1/d with d > 2, the local times only become well-defined once subtracted the 
divergent terms. This "renormalization" procedure motivates the study of a 
regularization. As a computationally simple regularization we discuss 

L Hf := J d 2 tS,(B H (t 2 ) - B„(i0), £ > 0, 



«B w ( f2 )-B ff ( fl )) : =(^ e- ' B "'-'J"'"" 2 



where 

Theorem 9 Let e > fre given. For all H G (0, 1) and a// dimensions d > 1 
t/ie functional L>h £ is a Hida distribution with kernel functions given by 



FH,e,2n(Ui, ...,U 2r , 

s B)7/ 4 ^T^rjn (MHl[, '* ])(Ui 



/or eac/i n = (m, ...,nd) G N d ; and Fh, £ ,yi = if at least one of the rij is an 
odd number. Moreover, if 2N(H — 1) + dH < I, then when e goes to zero 
the (truncated) functional Ljp e converges strongly in (S)* to the (truncated) 
local time L^p. 

Proof. As before, the first part of the proof follows from the Corollary |4] with 
respect to the Lebesgue measure on A. By the definition of the ^-transform, 
for all f G one finds 

S5 £ (B H (t 2 )-B H (t 1 ))(i) 

\ d 

1 \ g-i^z^yl/R^fW^ 1 ^,*,])^)! 2 



'2ir(e + t 2 - h 



2H\ 
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which is measurable. Hence, by Lemma El for all z G C and all f G <Sd(R) 
one has 

\S5 £ (B H (t 2 )-B H ( tl ))(zf)\< ( \ ==] e G ^\cXt'^\ 

\^2n(E+\t 2 -t l \ 2H ) ) 

with 2(^%~^-|2g) bounded on A and (e + \t 2 — ti\ 2H )~% integrable on A. 
By Corollary HI one may then conclude that Lu e G (S)* and, moreover, for 
every f = (/i, ...,/„) G S d ( 



SL H>E (f) = [ d 2 tS4(BH(t 2 )-B H (ti))(f) 

J A 



1 V f ^ 1 (_!)» 



2nJ J A (e+lh-t^")! ^ 2 n (e + \t 2 - h 

E n rr^n llf / ^/,(^)(M^i [tl)t2] )(x)) 

n l'"' ' n d j = l / 

nH hn d = n 



2inn 



v 7 ^; Ja 27 (e+^-tj 



12/nn+l 



2n, 

X) 



ni>— >»d j=l 
«lH n 



As in the proof of Proposition [HI it follows from the latter expression that the 
kernels F He il appearing in the chaos expansion L He = J2 n (: cu®" :,F Hen ) 
vanish if at least one of the in n = (ni, n d ) is an odd number, otherwise 
they are given by 

FH,e,2n(ui, ■■; «2n) 

To complete the proof amounts to check the convergence. For this purpose 
we shall use Corollary [31 Since 

SL%l(f)= [ d 2 *5'5f)(B if (t 2 )-B H (t 1 ))(f), 

J A 
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for all z G C and all f G S^R) we have 



SL^]{zi) < [ d 2 t \S8 £ (B H (t 2 ) - 'Bjj(tx))(zf)\ 

J A 



V2 



TVS 



V2 



Tit 



e ^ 



^tt| 2 | 2 ||f|| 2 



showing the boundedness condition. In addition, similar computations used 
to prove Theorem [7] yield for all {ti,t 2 ) G A 

\SS^(B B (t 2 )-BH(h))(t)\ 
1 \ d (C 2 



< 



< 



exPiV , -JL\t 2 - tl \ 2 - 2H \\f" 2 



' V ^t 2 -t 1 \ 2N ^- dH e^ f 



2 

H llf 1 1 2 



which allows the use of the Lebesgue dominated convergence theorem to infer 
the other condition needed for the application of Corollary [21 B 
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